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TECHNICAL NOTE D- 


RESONANCE- INERT STABILIZATION FOR SPACE STATIONS 

INTRODUCTION 


Space vehicles and space stations are rather flexible bodies whose 
resonances can destabilize attitude control systems, even though the control 
plant (i. e. , the flexible body) is passive. The controller consists of motion 
sensors (e.g. , platforms or rate gyros), signal processors, amplifiers, 
actuators, and control force devices (e.g., thrusters, gimbaled engines, 
torque motors, or control moment gyros). Regenerative undamping of the 
structure by the controller is possible and must be avoided through a good con- 
trol design. 

The same situation exists for the Skylab where the present attitude 
control system design depends also on a flexible body which must be properly 
described. The dynamic complexity of the body is overcome through an exten- 
sive analysis with progressive updating and a full-scale dynamic test of the 
short stack, which is used to verify and correct the modeling of the dynamically 
most-difficult Skylab section encompassing the Apollo Telescope Mount, 

Multiple Docking Adapter, and the Airlock Module as shown in Figure 1. The 
corrected model is then analytically combined with other parts, such as the 
S-IVB tank and the solar panels whose dynamic behavior is given by math 
models. All possible Skylab configurations (e.g, , undocked and docked) are 
thus obtained for designing the attitude control system. Obviously, present „ 
design methods do not take advantage of the flexible body's passivity. 

The method proposed herein departs from the conventional feedback 
loop concept and uses a controller with a response equivalent to a passive 
mechanical system consisting of springs, dashpots, and masses. The passivity 
of the controller yields a passive attitude control system which provides stiff- 
ness to the heading or the orientation of the vehicle and consequently satisfies 
the main control purpose, but also makes the system inert to resonances of the 
flexible body. The method was successfully used in 1966 in stabilizing a 
hydraulic support for dynamic tests with the Saturn- V space vehicle [ 1] at 
near free-flight conditions. 

The problem posed here is similar to the task of realizing certain filter 
characteristics by a passive electrical network, as it has been abundantly 
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Figure 1. Skylab cluster 










treated in the literature [ 2-14] for the past 40 years. The dynamic description 
of mechanical and electrical systems is analogous [ 15] , such that the electrical 
network theory is applicable to both cases, at least for one-dimensional prob- 
lems. Also, the analogy permits a discussion in mechanical terms, which are 
preferred here because control problems are thus directly interpreted. The 
simple feedback loop shown in Figure 2 serves as a brief explanation of the 
concept. The plant of Figure 2 is a passive mechanical system with a torsion 

spring, a , damping, a lf mass mo- 
ment of inertia, a 2 , and an additional 

'PLANT 

spring- mass system, S, where S can 
* also be described by a quotient. The 
controller consists of displacement 
and velocity sensors, their respective 
signal amplifications, b 0 and bj, and 



Figure 2. 


Feedback loop with lead-lag 
controller. 


the control force which lags with a 
time constant, r. If the system is 
constant (linear and time invariant) 
and initially at zero rest, then the 

response of the feedback of Figure 2 yields an equation in Laplace transform: 


M = 


</> 


/b n + sb, 

l 1 + ST 


+ a 0 


+ saj + s 2 a 2 




The factor on the right of equation (1) usually delivers the characteristic poly- 
nomial, whose roots determine the stability. Here, however, a reactive stabi- 
lization (discussed in the next three sections) will be applied which describes 
the structure of systems with two elementary connective operations, the addi- 
tion and the reduction 


a x 


a • b 

b = — 

a + b 


(2) 


For example, addition represents parallel arrangements of springs, where the 
total stiffness is the sum of the individual stiffnesses; reduction is a series 
connection of springs, where the total stiffness results from the reduction of 
the individual stiffnesses. These connections can also be used for dashpots and 
masses, where the latter is always connected to an inertial reference. The 
same connective operations are also found in electrical series and parallel 
circuits. 
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In equation (1), M is a moment and 4 > is an angular displacement; the 
objective is now to shape the response of the right side such that it behaves like 
a damped spring- mass system. Since a , a^, , and S represent a passive 

spring- mass arrangement, the system can be made stable, if the controller 
(b + sb )/ (1 + st) also behaves like a damped spring system. The main ques- 
tion now is how to structure the controller. With the help of a long division, we 
decompose first the quotient 


b 0 + sbi s (b 1 - b n r) 

1+ ST 0 i + ST 


( 3 ) 


into two additive parts and then with the help of the reduction into the structure 


b n + sbi 
1.+ sr 


= b n 


“ 


s(b 

l - b o t) 


( 4 ) 


Structure (4) is now completely decomposed such that mechanical elements can 
be recognized: spring b^, dashpot b^ - b^T, and again a spring (bj - b Q T)/T. 

The plus indicates that b is in parallel to the component contained in the brac- 
kets, where a reduction symbol indicates that b - b t and (b - br)/T are con- 

i o l o 

nected in series. However, structure (4) cannot consist of springs and dashpots 
unless the resolved parameters are all positive: 

4 

b , > 0, b - b t > 0, T > 0 . (5) 

v 1 0 


If inequalities (5) are satisfied, then structure (4) behaves like a damped 
spring system, which is not only stable and passive but also remains stable 
when connected to another spring- mass system (e.g. , S). A substitution in 
equation (1) yields 


M =. <p, 

b + 
J3 


(bi-b 0 T ) 

, x k_=blLl] 

1 T 

+ a +sa +s 2 a + S 

o l 2 

■ > (6) . 



- 

J 
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which with conditions (5) represents a damped spring-mass system as demon- 
strated in Figure 3. 



The system in Figure 2 is 
certainly stable, if conditions (5) are 
satisfied, which is intuitively clear 
from the mechanical equivalence in 
Figure 3. The mathematical proof is 
based on half-plane properties of the 
frequency response as defined and 
shown in the sections called reactive 
functions and reactive operations. 


If the system in Figure 2 is 

Figure 3. Mechanical equivalent of now treated the conventional way (e.g. , 
equation (6). after Hurwitz), then a characteristic 

polynomial is extracted from equation 
(1) where <p is the output. This means that S, which can also be a quotient, 
must be considered. Stating simply that S is passive will not suffice; all details 
of S must be known. If 


S= 0 


( 7 ) 


is assumed, then a third-order polynomial 


p(s) = a Q + b Q + s ^a t + a Q T + b^ + s 2 ^a 2 + a 1 s 3 a ? t 


( 8 ) 


is found, which results in the Hurwitz conditions 


a + a t 

a l 


a a T 


a + b 
o o 


a + a t + b 
i o i 


> 0, a 1 + a Q t + ^ >0 


( 9 ) 
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or the inequality 


b 0 a 2 T " a i (h + , a i T+a o 7 ") 

b > 

l a + a r 

2 i 


( 10 ) 


where all parameters are positive. 

The Hurwitz case gives a lower bound for the rate gain, , as the con- 
cept of mechanical equivalence which, according to inequalities (5), results in 
the simple condition 


b > b r . (11) 

10 


A comparison of inequalities (10) and (11) shows that mechanical equiva- 
lence is more conservative than the conventional approach. But mechanical 
equivalence accepts any passive plant with numerous and unknown resonances, 
while the conventional approach requires an exact specification of the plant 
(e.g. , S = 0) as in this example. If the plant has no damping (a^ = 0) then 

Hurwitz condition (10) and inequality (11) become equally restrictive. 

Mechanical equivalence is based on the following reactive stabilization 
principle. 

Definition 1: A system, connected by addition, is reactively stable when 
the system is stable in the decomposed and the assembled state. 

Reactive stabilization is similar to the absolute stability of electrical 
networks [ 16] which are stable regardless of terminating impedances. Here, 
the elements or building blocks are reactive functions as defined in the next 
section, structural connections are reactive operations as given in the second 
section, and n-dimensional presentations are treated with reactive matrices in 
the third section. 

Stability analysis on subsystems (here, the controller) is thus possible 

and plant changes caused by different passive configurations-do .not-affect - 

Stability; this apparently is an advantage in docking and deployment situations 
of space stations. 
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REACTIVE FUNCTIONS 


Reactive stabilization is applicable to constant systems and represents 
systems or subsystems by mappings of the right half-plane's contour with 
Nyquist-type plots [ 17] relative to the origin. Reactive connections, [i. e. , 
additions and reductions (equation (2)] are important because they propagate 
half-plane properties of mappings of the contour domain depicted in Figure 4. 

The contour about the first 
quadrant can also be used to enclose 

the fourth quadrant by taking the 
complex conjugate. Any poles and 
zeros, enclosed in the right half- 
plane and the contour of the first 
quadrant, are thus related to the map- 
ping^ s encirclements of the origin in 
multiples of 180 deg (N ). From the 

7T 

residue theorem [ 18, 19] follows 

N = Z - P, (12) 

7 r 


where Z and P are the zeros and 
poles, respectively, within the contour of Figure 4 and the contour’s complex 
conjugate. Figure 4 describes one domain for the contour mapping; however, 
three mapping ranges were selected according to Figures 5 through 7 and the 
definitions of three different reactive functions: positive imaginary, positive 
real, and negative imaginary. 

The domain and the ranges of Figures 4 through 7 are defined as 
follows. . 

Definition 2: The domain of contour mappings by reactive functions is 
the first quadrant with a closed bound on the positive real axis and open bounds 
at the origin, the point of infinity, and the positive imaginary axis (Fig. 4). 

Definition 3: The range of contour mappings of positive imaginary 
(denoted + j) functions is the upper half-plane with a closed bound on the posi- 
tive real axis and open bounds at the origin, thepoint of infinity, and the nega- 
tive real axis (Fig. 5). 


1 +j 



Figure 4. Contour of the first quadrant 
with p > 0, lim p= 0, e> 0, lim e = 0, 
co> R>0, and lim R = ® . 
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Figure 5. Range of contour mappings for positive imaginary functions 

with example s 2 . 


OPEN BOUND CLOSED BOUND 



Figure 6_. Range of contour, mappings for- negative imaginary 

-2 

functions with example s 
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Definition 4: The range of con- 
tour mappings of positive real (denoted + 
r) functions is the right half-plane with 
open bounds at the origin, the point of 
infinity, and the whole imaginary axis 
(Fig. 7). 

Definition 5: The range of con- 
tour mappings of negative imaginary 
(denoted - j) functions is the lower 
half-plane with a closed bound on the 
positive real axis and open bounds at 
the origin, the point of infinity, and the 
negative real axis (Fig. 6). The point 
of infinity considered here is one pole 
on the Riemann sphere [18, 19] . 

Definition 6: Functions are 
reactive iff 1 they have no poles within 
the contour of definition 2 and the con- 
tour’s complex conjugate, iff they map 
the positive real axis on itself, and iff they map the contour of definition 2 
within one of the ranges of definitions 3 through 5. 

Theorem 1: All reactive functions imply no zeros within the contour of 
definition 2 and the contour's complex conjugate. 

Theorem 1 is proven in denying the implication and assuming that at 
least one zero (Z = 1) exists. But the contour mappings of definitions 3 through 
5 do not encircle the origin (N =0), reactive functions of definition 6 do not 

7T 

contain any poles within the contour of definition 2 (P = 0), and the hypothesis 
of reactive functions yields with equation (12) a Z = 0. Therefore, the denial 
of the implication (Z = 1) leads to a contradiction and theorem 1 is proven. 

Theorem 2: All (+ j)-functions imply for s = 0 and s = « a limit cs n 
(c > 0), where n e { 0, 1, 2} . 

Theorem 3: All (+ r)-functions imply for s = 0 and s = ® a limit cs 11 
(c > 0), where n e (-1, 0, 1} . 


1+i 



Figure 7. Range of contour mappings 
of positive real functions with example 
s + 1 /s . 


i . If and only if. 
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1 : Theorem~4: All (- j) -junctions imply for s = 0 arid s = 00 a limit cs 11 

(c > 0) , where n e (-2, -i, 0}. 

Theorems 2 through 4 are proven by denying the implications which 
means in the order of the theorems that n /{0, 1, 2};' n /{-i, 0, l}, and 
n -i, O}. A substitution of the boundaries p > 0, lim p = 0, e> 0, 

lim e=0, 00 > R > 0, and lim R = 00 into s' 1 leads to contradictions with reac- 
tive, functions as given by definitions 2 through 5. The open bounds of Figure 4. 

are represented for the origin by a small quarter circle with s = pe^ 5 (p > 0; 
lim p = 0, 0^5<7 t/ 2), for the point of infinity by a large quarter circle with 

s ,= Re^ (°° > R > 0, lim R = <*>, 0 ^ 6 <?/%)> and for the positive imaginary axis 
by a line with s = e + jo> ( e > 0, lim e = 0, pSusR)., A, substitution of these 

boundaries into s 11 proves theorems 7 through 9. 

Theorem 5: All '(+ j) -functions q(s) . excluding constants, imply a 
mapping of the line s = e+ jw (e > 0, lim e = 0, p £ < R) such that 

Im [q( e + jw) ] > 0. 

Theorem 6: All (+ r) -functions q (s) imply a mapping of the line 
s =■ e + jw(e > o, lim e = 0, p ^ co < R) such that Re [q(e + jo;)] > o. 

Theorem 7: All (-j)-functions q(s), excluding constants, imply a map- 
ping of the line s = e + jcn (e > 0, lim e = 0, p < cn <R) such that 
Im [ q(e + jco )] < 0. 

Theorems 5 through 7 are proved by denying again the implications 
which produce in the order of the theorems Im [ q(e + jcu ] <0, Re [ q(e + jo> )] <D, 
and Im [q(e + jcu )] >0. Such mappings would lie outside the range specified 
in definitions 3 through 5, thus contradicting the hypothesis of reactive func- 
tions, and theorems 5 through 7 are proven. 

Theorem 8: A function q(s) is a reactive iff q(s) is stable in the sense of 
definition 6, iff all coefficients are positive real, and iff at least one of the 
theorem-pairs 2 and 5, 3 and 6, or 4 and 7 holds. 

The proof begins with the "only if" part which equates reactive functions 
with definitions 2 through 6. The latter produces all theorems (2 through 7) 



used in theorem 8. Definition 6 implies (asymptotic or semi-) stability and 
positive real coefficients. Definitions 2 through 5 imply that the limit cases 
(s = 0, s = 00 ) and the mapping of the open bound on the positive imaginary axis 
must be in the same half-plane, which results in the pairing of theorems 2 
and 5, 3 and 6, and 4 and 7. Definitions 3 through 5 give three choices, imply- 
ing at least one of the three pairs must be valid. 

The proof continues with the "if" part. Stability and positive coefficients 
imply definition 6. The pairing of theorems 2 through 7 implies that the ‘limit 
cases of theorems 2 through 4 must lie in the same half -plane as the mapping 
of the open bound on the positive imaginary axis (theorems 5 through 7) . 
Therefore, the contour of definition 2 is mapped at least within one of the ranges 
of definitions 3 through 5. Thus, definitions 3 through 6 are met which implies 
that the function q(s) is a reactive function. Both proof parts show that 
theorems gives necessary and sufficient conditions for reactive functions. 


The stability of reactive functions is found directly from the denomina- 
tor's polynomial which is also related to the following mapping property. . 

Theorem 9: Reactive (+ j)-, (+ r)-, and (- j) -functions imply respec- 
tively that the minimum of all imaginary parts, the minimum of all real parts, 
and the maximum of all imaginary parts lie only on the contour of definition 2 . 


The proof of theorem 9 follows from the maximum principle of harmonic 
functions [19] which relates the analytic property of these functions to the 
maximum at the boundary of a simply connected s-domain. For example, in 
the case of a (+ j) -function q(s) = u(s) +jv(s) [v(s) ^0], we use the auxiliary 
form 


w(s) = e 


jq(s) 


_ -v (s) ju (s ) 
- e e 


(13) 


and find from the differential 


w' (s) = jq' (s)e^ S ^ 


(14) 


that function w (s) is analytic if the function q (s) is analytic and satisfies the 
Cauchy-Riemann differential equations [18, 19]. Since reactive functions q(s) 
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are analytic, w(s) must be analytic. If we deny now the implication for (+ j) - 
functions and state that the imaginary part's minimum v(a) = m (m > 0) lies 
at a point "a" inside the contour of definition 2, then the analyticity yields with 
Gauss' mean -value theorem the inequality 



x 

2 ?r 


/V T( a + r eW d<f> 
0 


(15) 


Since -m was assumed to be a maximum, there must be at least one angular 
section A <p - (f) 2 - where 


-m ~v(a + rej<£) -m-e 
e > e = e 


e > 0, and <p i < (j> < <p 2 


(16) 


holds. Consequently, inequality (15) results in 


-m _ -m 
e ^ e 


1 - 


A 4> 
2n 




< e 


-m 


(17) 


which is obviously a contradiction. Therefore, the maximum cannot be at a 
point "a" interior to the contour of definition 2, and theorem 9 is proven for 
(+ j ) -functions . The proofs for (+ r)- and (- j) -functions are similar, if the 

respective functions w (s) = e = e e and w (s) = e 

v (s ) — ju ( s ) 


= e 


are used. 


Theorem 10: A function q(s) is reactive iff real values s = r are mapped 
as real values q(r) and iff one of the ranges of definitions 3 through 5 is 
occupied by the mappings of all the contours within the contour of definition 2. 

The proof begins with the "only if" part that a function is reactive which, 
according to definition 6, implies that real values s = r are mapped as real 
values q(r) . Reactive functions also map the contour of definition 2 within the 

ranges specified by one of definitions 3 through 5. Definition 6 and theorem 1 

show that the con to_ur_of_definition-2-contains no~pole~s~ahd~zeros: therefore, 
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mappings of any contours interior to definition 2 will not encircle the origin 
and with theorem 9 will stay within one of the ranges of definitions 3 through 5. 
Therefore, all implications of the "only if" part are met. 

The second half of the proof continues with the "if" part of mapping the 
positive real axis onto itself; this satisfies the real part statement of definition6. 
The "if" part of the ranges satisfies also definitions 3 through 5 for all contours 
within definition 2, which implies that the function does not have any poles and 
zeros within the contour of definition 2 . All conditions of definition 6 are met 
and the function is reactive. Both halves of the proof are thus completed and 
theorem 10 gives necessary and sufficient conditions. 

Well known in network theory [2-14] are positive real functions which 
are the (+ r) -functions of definition 4 when definition 6 is also satisfied. 

Definition 7: A function q(s) is positive real iff real values s = r result 
in Im [q(r)]= 0 and iff s e [s; Re (s) > 0] produces Re [q(s) ] >0. 

Definition 7 and theorem 9 result in the properties of definitions 2, 4, 
and 6. Since definition 7 and theorem 9 are not easily applied, the contour 
approach of definitions 2 through 6 is preferred here. Note also that (+ r)- 
functions result in (+ j) - or (- j) -functions when multiplied by s or when 
divided by s, respectively. 

Definitions 3 through 5 and Figures 5 through 7 represent three different 
types of reactive functions, depending on the half -plane property of a transfer 
function or a quotient. For example, a constant element is a (+ j)-, (+r)-, and 
(- j) -function; an element s is a (+ j) - and a (+ r) -function; an element s -1 is 
a (+ r)- and a (- j) -function; an element s +2 is a (+ j) -function; an element 
s + s -1 is a (+ r) -function; and an element s -2 is a (- j) -function. These half- 
plane properties are preserved when the elements are properly connected and 
thus permit the construction of stable systems because the resulting reactive 
functions have no poles and no zeros within the contour of definition 2 accord- 
ing to definition 6 and theorem i . 


REACTIVE OPERATIONS 


Reactive stabilization is a building -block approach, where the elements 
belong to one reactive function set, specified by definitions 3 through 5, and 
where the structure is given by reactive operations. 
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Definition 8: Operations producing reactive functions from functions of 
the same reactive type are denoted as reactive operations. 

If an operation of definition 8 connects two (+ j) -functions, a (+j) -function 
will result. The example applies also to (+ r)- and (- j) -functions; but mixing 
between the three function sets is not allowed here and will not be discussed, 
even through mixing is advantageous in certain cases [ i ] . 

Theorem ii: The addition is a reactive operation. 

Theorem 12: The reduction, defined in equation (2), is a reactive 
operation. 

improving theorem 11 we observe that the mappings of the contour from 1 
definition 2 are contours which are described by a complex vector. Therefore, 
the addition of reactive functions is a vector addition where the resultant is 
given by the addition of the vector's components. From the half -plane proper- 
ties reflected in theorems 2 through 7 and definitions 3 through 5, we see that the 
vectors are either positive real or complex. When complex, the vectors have 
at least one (either positive or negative) definite component. The addition of 
positive real parts results always in positive real parts and the addition of 
positive (negative) imaginary parts will always produce positive (negative) 
imaginary parts. Therefore, the addition preserves the half -plane property of 
reactive functions as long as the function types are not mixed. Also, the addi- 
tion preserves stability because only stable elements are added. Thus, all 
conditions for reactive functions (definitions 2 through 6) are met, definition 8 
is satisfied, and theorem 11 is proven. 

For the proof of theorem 12 we consider first the reductive vector 
presentation 


z = a x jb = a — 


a + b 


7 + jb t , , 

1 a* + b 


2 ~ U + JV 


(18) 


where a, b, u, and v are real numbers. Equation (18) shows that the reduc- 
tive and the additive components are of the same polarity, which means that 
the positive or negative definitiveness of one component, as given by a reactive 
function set, will not be changed by the choice of the vector presentation. For 
reductive conne ctions it is advarit ageou s_to _u seethe .r.educ tive -presentation : — 
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z = Z! x z 2 = a* x jbi x a 2 x jb 2 = (a t x a 2 ) x j (b t x b 2 ) = a x jb , (19) 


because the resultant components 


a = ai x a 2 and b = b A x b 2 


( 20 ) 


are simply determined. Equation (2) shows that the reduction does not change 
the polarity when real numbers with the same sign are reduced; e.g. , if b^ 
and b 2 are both positive (negative) definite, fheir reduction will also be posi- 
tive (negative) definite. Thus, the half-plane characteristic of reactive func- 
tions is also preserved by this operation. . Further, we find that the reduction 
of two reactive quotients - 


Ni(s) N,(s) ^(s) - No (s) 

Dj(s) D 2 (s) D t (s)N 2 (s) + D 2 (s)N 1 (s) 

produces a stable common numerator N i N 2 with only stable zeros (Z = 0) , 
since the individual numerators N t and N 2 have only stable roots (theorem 1) . 
The stability of the zeros (Z = 0) and the half-plane property (N = 0) causes 

7T 

only stable poles in equation (21) because equation (12) results in P = Z = 0. ’ 
All requirements of reactive functions (definitions 2 through 6) are thus met, 
definition 8 is satisfied, and theorem 12 is proven. 

Analogous to the subtraction we introduce now the negative reduction or 
the subdue ti on 


a \b 


ab 

b - a 


( 22 ) 


which together with the subtraction is not a reactive operation. Between the 
addition and the reduction exists a duality which is listed in Table 1 . The 
numbers a, b, and c are complex and are unequal to zero or infinity, where 
infinity represents the point of infinity on the Riemann sphere [18, 19 ]. 

Addition and reduction form two Abelian groups which are coupled through the 
unification rule. Both groups are not distributive among each other, but are 
distributive with multiplication and form two rings [20] which are again coupled 
through the unification. 
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TABLE i. DUALITY OF ADDITION AND REDUCTION 


Associativity: 

a + (b + c) = (a+c) + c 

a x (bxc) = (axb) xc 

Distributivity: 

a (b + c ) = ab + ac 

a (b x c ) = ab x ac 

Commutativity: 

a + b = b + a 

a x b = b x a 

Repetition: 

a + a + a = 3a 

a x a x a = a/3 

Perpetuation: 

a + a + a+ ... = °° 

a x a x a x . . , = 0 

Identity: 

a + 0 = a 

a x °°= a 

Cancellation: 

a - a = 0 

a \ a = 00 

Unification: 

(a + a) x (a + a) = a 

(a x a) + (a x a) = a 

Inversion: 

(a + b) -1 - a -1 xb' 1 

(a x b) -1 - a -1 + b _1 


Multiplication cannot be identified as a reactive operation, since it can 
rotate a complex vector outside the range of a specified half-plane. In special 
cases, multiplication produces reactive functions, but generally this is not 
true. However, the inversion always preserves reactive functions. 

Theorem 13: Inverse reactive functions are reactive functions. 

Theorem 14: Inversions are reactive (+ r) -operations. 

Inversions produce vectors pointing in the direction of the complex 
conjugate and thus do not change the real part, but they reverse the polarity of 
the imaginary part. The half -plane properties are preserved for (+ r) - . 
functions, but (+ j) -functions exchange the half -planes with (- j ) -functions and 
conversely. Poles and zeros are only rotated, but are stable for reactive 
functions. Thus, definitions 3 through 6 are satisfied and theorem 13 is proven. 
For (+ r) -functions, definition 8 also holds and theorem 14 is proven. 

The synthetic approach to stability was first published in 1966 [1] with 

a different nota tion (e.g., reduction was then called inver se addition ),._In 

1967 [21], Lewis suggested a complementary algebra for describing electrical 
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switching circuits, where he denotes a "complementary addition" or "sulp" for 
the operation which is called reduction here. The algebraic features of systems 
are obviously present and can be a valuable tool if properly defined as in the 
case of the addition and reduction which propagate stability. 


REACTIVE MATRICES 

In multivariable feedback cases [22] inputs and outputs are related by a 

transfer matrix A(s) whose elements a., (s) are transfer functions or quotients. 

lk 

Also, an attitude control system of a space station belongs to this category, 
because of the three-dimensional feedback from the structure as shown for the 
Skylab in the next section. This section brings the theory of reactive stability 
for n-dimensional cases based on quadratic forms y*Ay (asterisk means complex 
conjugate transpose). 

Definition 9: A transfer matrix A(s) is stable, when all elements a.^s) 
are (asymptotic or semi-) stable transfer functions. 

Definition 10: Matrix A(s) is complex definite, when the quadratic form 
q(s) = y* Ay has for all vectors y ^ 0 or 00 at least one complex component 
which is either positive or negative definite. 

Definition 11: Matrix A(s) is reactive, when A is complex definite and 
nonsingular, and when the quadratic form q(s) = y* Ay is a reactive function 
for all vectors y * 0 or 00 . 

Theorem 15: Reactive matrices are stable. 

The proof uses definition 11, where the quadratic form y*Ay is a 
reactive function for all vectors y^ 0 or °°. The quadratic form combines 
linearly all matrix elements -a (s) and therefore draws all its poles from 

these elements, except for possible pole-zero cancellations. General cancella- 
tions involve a factor k common to all matrix elements a., = k a., , such that k 

lk lk 

is still a part of each element; special cancellations for certain y-vectors can- 
not hide any poles because other y-vectors will show them. Therefore, the 
quadratic form y* Ay for all vectors y ^ 0 or 00 represents all poles of the ele- 
ments a which must be stable by definition 6, since y* Ay is a reactive func- 
tion. With definition 9, stable elements lead to a stable matrix, and theorem 
15 is proven. 
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Theorem ,16: .Nonsingular complex definite matrices have complex defi- 
nite inverses. 

The proof uses the transformation . 


z = Ay 


(23) 


and the quadratic form 


q* (s) = y* A* y = z* A _1 z 


(24) 


By hypothesis, matrix A was complex definite. Since the complex conju- 
gate of the quadratic form q* = (y* Ay)* = y* A* y only changes the sign of 
the imaginary part, matrix A* must be also complex definite. Since all vectors 
y / 0 or 00 are involved (definition 10) and matrix A is nonsingular by hypothesis, 
transformation (23) produces all possible vectors z f 0 or °°. Equation (24) is 
therefore valid for all possible vectors y and z , and the quadratic form 
q* = z* A -1 z must be complex definite for all vectors z ¥■ 0 or thus j proving 
theorem 16. 

Theorem 17: Inverse reactive matrices are reactive. 

The proof begins with the form y* A* y of equation (24) which is a 
reactive function for a mapping of the complex conjugate contour of definition 
2. This means that y * A* y has the half -plane and the stability properties of 
reactive functions. Theorem 10 is therefore applicable even though complex 
conjugate contours are used here .• All contours within the complex conjugate 
contour of definition 2 must result in one and the same half -plane property for 
all vectors y^O or “. Since transformation (23) is not singular and equation 
(24) holds, the quadratic form z* A -1 z must maintain the same half-plane 
properties as y* A*y for any vector z ¥ 0 or 00 . Additionally, all real values 
s = r are mapped as real values q* (r), thus satisfying all conditions of theorem 10 
and the form z* A -1 z is a reactive function for any vectors z ^ 0 or °°. There- 
fore, matrix A -1 is by definition 10 complex definite, by. assumption nonsingu- 
lar, and consequently reactive, by definition 11. This proves theorem 17. 

Note that by theorem 15, matrix A -1 is also stable, which is important in 

proving overall stability for n-dimensional cases, especially when c onnect ions 

are .used-as-emphasized by-the next statementsr ~ “ 
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Definition 12: The reduction of matrices A and B is defined as 
A x B = (A -1 + B' 1 ) -1 . 

Definition 13: Operations producing reactive matrices from matrices 
of the same reactive type are denoted as reactive matrix operations. 

Theorem 18: Addition is a reactive matrix operation. 

Theorem 19: Reduction, according to definition 12, is a reactive matrix 
operation. - 

The proof of theorem 18 uses the quadratic form 


y* (A, + B) y = y* Ay + y* By , 


(25) 


where A and B are reactive matrices. The right of equation (25) is a sum of 
reactive functions which again is a reactive function. This is valid for all 
vectors y ^ 0 or °°. Equation (25) equates all properties of reactive matrices 
with the quadratic form of matrix A + B which therefore must be reactive, and 
theorem 18 is proven. 

The proof of theorem 19 employs the quadratic form 


z* (A -1 + B *)z — z* A 1 z + z*B~*z 


(26) 


where A and B are reactive matrices. By theorem 17, A -1 and B -1 are also 
reactive matrices and by theorem 18, A -1 + B -1 is also a reactive matrix. 
Definition 12 leads to the quadratic form 


y* (A xB)y= y* (A -1 + B -1 )""V 


(27) 


where A -1 + B _1 is a reactive matrix whose inverse is by theorem 17 reactive 
again. Therefore, the reduction of matrices of the same reactive kind pre- 
serve the reactive property, definition 13 is satisfied, and theorem 19 is 
proven. Note that the reduction represents a difficult matrix operations, but 
permits a rather easy stability assessment if matrices A and B are reactive. 
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ATTITUDE CONTROL LOOP 


The following sections demonstrate that reactively stabilized attitude 
control systems are feasible for space stations. The Skylab's control moment 
gyro (CMG) loop [23-27] has been chosen to show a possible solution with rate 
gyros and angular accelerometers as attitude sensors. The sensors are 
mounted near the CMGs on the racks of the Apollo Telescope Mount (ATM) as 
depicted in Figure 8. The discussion is limited to the CMG loop, even though 
the principle can be extended to include the experiment pointing control system 
(EPCS) and the thruster attitude control system (TACS). The CMG control 
loop is shown in Figure 9 by a vector signal flow diagram. 


The commanded attitude is given by the vector cp, C is a control gain 
matrix, G is the CMG response matrix, M is the moment vector applied to 
the Skylab structure, P is the Skylab's structural response matrix, y is the 
attitude vector measured by three rate gyros and three angular accelerometers, 
and Q is the sensor response matrix feeding back to the control gain matrix G. 
Because the Skylab rotates slowly, coriolis and centrifugal forces are neglected 
and the dynamic behavior is expressed by the linearized vector form 


7 = PM , (28) 

The plant response is given in Laplace transform by the transfer matrix P 
which is the sum 



+ K 


(29) 


of the moment of inertia matrix 



’ll 

0 0 


V 

J = 

0 

I 2 0 

- 

h 


0 

o Is. 


J 3 _ 


[Ji J 2 J 3 ] 


(30) 
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Figure 9. Vector signal flow diagram of a CMG 
- attitude control system. 

and the resonance matrix 


K = Y'fi -1 Y 


(31) 


Matrix K is the product of the transpose of the mode matrix 



(' means transposed) and the modal resonance matrix 


(32) 
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0 


0 


mj (s 2 + s2 w[) 0 

0 m 2 (s 2 + s 2 ^ 2 w 2 + °4) " 0 ' . . . 0 


0 0 m 3 (s 2 + s2£, 3 o^ + o;|) 


0 


12 = . 


• (33) 




0 0 0 m (s 2 + s2£ w + cJ ) 

n n n n 


The inertial matrix, equation (30) , is always the sum of a diagonal and a diadic 
matrix. The resultant elements are moments and produce moments of inertia. 
The mode matrix elements, equation (32), are modal slopes at the area of the 
CMGs and the nearby located attitude sensors (rate gyros and angular accelero- 
meters) . The first index of the elements represents body -fixed coordinates i, 

2, or 3, and the second index is related to the modal number. The modal 
resonance matrix has only diagonal elements which represent resonators with 
generalized masses m., damping factors and resonant frequencies ^. 


QUADRATIC FORMS 

The passivity of the Skylab structure follows from the quadratic form 

. M* y = M* PM = ~ M* J -1 M + M* Y’ fi _1 YM = q t , (34) 

s 

where M* is the conjugate transpose of the vector M. Since the inertial matrix 
J is always positive definite, we find. with M = Jz that 


M* J _1 M = z* Jz (35) 
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and therefore the form M* J -1 M is positive definite. Since (YM)* = M* Y', 
the form (YM)* fi _1 (YM) results in a positive sum of second-order systems. 
Therefore, matrix P is complex definite (definition 10) and has a quadratic 
form which is a reactive (-j) -function (Fig. 6, definitions 2, 5, and 6) . 
Consequently, matrix P is reactive as specified in definition ii. Since the 
plant P is a reactive (- j) -matrix, the controller must also be a reactive (- j)- 
matrix. The system in Figure 8 now yields the vector equation 


<p = (QP + C -i G _1 ) M 


(36) 


which is changed to 


Q -1 <(> = (P + Q -1 C _1 G -1 )M 


(37) 


A new quadratic form 


q(s) = M* (P + Q -1 C -1 G -1 )M = qi + qa 


(38) 


results where q t is already given by equation (34) and q2 is defined by 

q2 = M* Q _l C" 1 G" 1 M (39) 

Equation (39) is the quadratic form of the controller whose reactive stability 
property must be further investigated to ensure that q 2 becomes a reactive 
(~ j) -function. This principle permits a separation of the controller from the 
plant for a stability analysis similar to the example of the introduction 
(Figs. 2 and 3) . , 


CONTROLLER RESPONSE 


The response of the CMGs 2 is given by a diagonal matrix 


1 + s2f, /oj 0 + s 2 / co 2 

2. Data contributed by H. L. Shelton, S&E-ASTR-SD 


( 40 ) 


24 



where I is a 3 by 3 identity matrix. The angular sensor response is assumed 
to be negligible as expressed by the matrix 


Q = I 


(41) 


The control gain matrix is 


C = 


c 01 + sc u + s^c 2 i 0 0 

0 c 02 + SC 12 + s 2 c 22 0 


0 


0 


c 03 + sc 13 + sm 23 


(42) 


where e 0 j is the integrated rate gyro gain, Cjj is the rate gyro gain, and c 2 j is 
the angular accelerometer gain; therefore, the total block (GCQ) -1 has the. gain 


Q^C^G' 1 = 


921 0 0 


0 922 0 


0 0 C123 


1 + s2£, / 0J n + S 2 / 

<l21 C 01 + SCu + S 2 C 2 1 

_ 1 + s2E, / + s 2 / 

^ c 02 + sc 12 + s 2 c 22 

_ i + s2t, / + s 2 / a; 2 

^ c 03 + sc 13 + s z c 23 


(43) 


DECOMPOSITION OF CONTROLLER RESPONSE 

The controller is represented by the quadratic form of equation (39) 
which with equation (43) produces the sum 


92 = | Mi | 2 q 2 ! + | M 2 | 2 q 22 + |m 3 | 2 q 23 (44) 

where all quotients q 21 have only positive factors M. 2 , Therefore, q 2 will be a 

(- j) -function if all q 2 j are reactive (- j) -functions like q l5 equation (34) . The 
shaping is found by a decomposition into reactive (- j) -elements which are 
reactively connected (definition 8). First, we separate from q 2 j a constant 
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Pl c 2i a) 0 


(45) 


and obtain a remainder 


r - n n - 1 " CniPl + S (2 ^ /W 0 ~ CiiPi) _ 

r t - 021 Pl „ . + r 2 X r 3 

c 0l + SCii + S C 2l 


(46) 


which is split by a reduction into two components 


r 2 


= 1 - CpiPi + s (24/<a 0 - c liPl ) = 


s 2 c 


r 4 + r 5 


2i 


(47) 


and 


r3 = - 1 -£flLPj L±li 2 i A^l L JL£liP l > = r e + r 7 
c oi + scjx 1 


(48) 


From r 2 and i' 3 , we obtain the new components 


„ _ 1 “ c oiPl - P2 
4 ~~s z c^i~~ s 7 


(49) 


A 


gjAfr ~ onPi = £3^ 


SC 2 i 


(50) 


- 2^/^q -c liPl ^ 

r 6 . P4 

C 11 - 


(51) 


and 


~ „ - 1 c oiPi “ c oi r 6 — s. „ 

r? 3 - r e — + — r 8 X r 9 

Cfli 1- SCjx 


(52) 
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The remainder r 7 is decomposed further by a reduction which produces the 
elements 


r 8 = 


1 - 


c oiPl ~ c oi r 6 _ 


c oi 


Ps 


( 53 ) 


and 


r = 1 c oiPl c ni r fi _ Ps 
9 SCii s 


(54) 


This completes the decomposition which yields the following structure 


<l2i = Pi + 





(55) 


where the elements have factors s” with n e (-2, -1, 0} as given by equations 
(45), (49), (50), (51), (53) , and (54) . If all parameters p i are positive, 

then q 2 i are (- j) -functions and q 2 is a (- j) -function, because reactive opera- 
tions connect (- j) -function elements. 


REACTIVE STABILITY CONDITIONS 


Decomposition (55) results in the same conditions as required to make 
the imaginary part of q^jw) negative. Equation (43) and theorem 7 yield for 
oo > 0 the condition 


Im[q 2 i(jw) ] = w 


f-Mcpi-^ca) - c lt (i -;j) 

(c oi - w*c 2i ) 2 + Jofi 


< 0 


(56) 


Additionally, if all coefficients (c 0 i, cjj, c 2 i, f,, and tx> 0 ) are positive real, 
then the limits for s = 0 and s = 00 will be positive constants, and the poles 
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and zeros will also be stable. Theorems 4, 7, and 8 are satisfied and all q^ 
are reactive (- j)-functions. Inequality (56) results in the negative definite 
polynomial 




u> 


0 


c oi ~ c li + w 2 





(57) 


which yields the inequalities 


c 2i a, o 


utf > 


Ci i OJ, 


2t, 


> Cgi > 0, 


i e{l, 2, 3} 


(58) 


The same conditions are derived from decomposition (55), which is therefore 
a necessary and sufficient decomposition. Often, decompositions yield only 
sufficient- (conservative) conditions which are unnecessarily restrictive and 
therefore must be avoided; e.g., by a proper structural pattern. Since the 
decomposition is rather difficult, an evaluation by complex component parts is 
recommended. 

Inequalities (58) produce for Skylab's CMG control system co 0 = 34 and 
t, = 0.5 the following reactive stability condition: 


c 2 i 1160 > Cjj 34 > c 0 j > 0 


(59) 


The coefficients c 0 i and c^ are determined from the rigid-body resonances 
which are found by an approximation for small s values. Equations (28), (29) , 
and (36) yield the rigid -body response 


C (p = (C + s 2 J) y 
with the control matrix 


(60) 


c 0 i + sc tl 0 0 

C = 0 -c 02 -* sc t2 

0 0 c„ 3 + sc 13 


( 61 ) 
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and the inertial matrix of equation (30) . The Skylab's inertial matrix [25] in 
kgm 2 is 


101.6 -2.155 -38.4 


J • 10" 4 = 


-2.155 


-38.4 


615 -3.45 


-3.45 


605 


125.6 

0 

0 


~ 4. 9 ~ 

0 

615.2 

0 

- 

0.44 

0 

0 

666.5 


7. 8 


[4.9 


0.44 


7.8] 


(62) 


Equation (60) has the form 


/ 

n " 


r- - 

V 

C <t> =/ 

c 0l + SCi! + S^Ij 0 0 

0 c 02 + sc 12 + s 2 I 2 0 

-s 2 

Jl 

J 2 

|^1 J 2 ^ ) 

( 

0 0 Cq 3 + sc 13 + s 2 I 3 


_J 3 _ 

) 


(63) 


where the diagonal elements are assumed to be proportionally designed: 


c 0 i + scti + s 2 Ii = -± (c 0 + scj + s 2 I 0 ) , i e {1, 2, 3} 

i o 

such that equation (63) takes the simple form 



'll 

0 

0 ' 


Ji' 


c 0 + sc t + S 2 I 0 

0 

I 2 

0 

_. q 2 

T 

Tt t t! 

Io 

b 

J 2 

1_ J 1 "2 J 3J' 


0 

0 

1 3. 


_ J 3_ 



(64) 


(65) 
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The inverse of equation (65) is now 


( 

J1A1 

7 =(I + 

<^2/^2 

f 

J3/I3 


[ J 1 J 2 J3 I S 2 I 0 
c 0 + scj + sT 0 (1 - n) 


c 0 + SC! 
c 0 + SCj + s 2 I 0 


0 


with the abbreviation 


(66) 


II 



+ 


t2 t2 

" 2 . + 

I 2 J 3 


(67) 


Equation ( 66 ) relates the commanded orientation <p with the space station's 
angular position y and shows two undamped resonances 



and 

= JiTu^i7 (68) 

which are rather close. For example, Skylab has a II = 0.285 and the 
resonances differ only by approximately 14 percent. 

The following approximate parameters are assumed for Skylab: 
c 01 = 16 000 Nm, Cq 2 = 79 000 Nm, c 0 3 = 85 000 Nm, and c t i = 10 c 0 i Nms. Now 
inequalities (59) yield the reactive stability condition 

c 21 1.16 > 5400 > 16 > 0 

c 22 1.16 > 27 000 > 79 > 0 — (69') _ 

c 2 3 1.16 > 29 000 > 85 > 0 


30 



The Skylab case satisfies all conditions except for the c 2 i's which are zero, 
since no angular accelerometers are used - Skylab has been conventionally 
stabilized by attenuating sufficiently the controller's response at structural 
resonances. The lowest rigid -body resonance is approximately 

f i0 = 0.018 Hz, Cm = 0.6 . (70) 


The present control system of Skylab uses the conventional approach 
with rather low control gains when compared to launch vehicles. However, 
these gains give the Skylab sufficient attitude stiffness and damping against 
disturbances; e.g. , crew motion. The control gains are changed depending on 
the docking of the Apollo command and service module. Since reactively 
stabilized systems are insensitive to docking configuration, the new concept 
appears to be ideally suited for such applications; however, scheduling problems 
make any new design changes prohibitive, especially when the conventional 
approach has already produced a workable, attitude control system. 


CONCLUSION 

This report demonstrates that system structures can make controllers 
inert to plant resonances if the principle of reactive stabilization is applied. 
Overall stability can thus be predicted by analyzing only the controller without 
knowing the resonances of the passive plant. This is especially advantageous 
for space stations (Fig. 10), where configurations vary widely with deployment 
and docking and where conventional methods may be impractical because of the 
analysis and bookkeeping problems associated with flexible body resonances. 

Today's sophisticated computer programs have only partially succeeded 
in math -modeling dynamic behavior of complex systems, and full-scale dynamic 
tests are still used for model verification. Also, structural damping is difficult 
to predict and must usually be determined by tests. Small damping factors 
cause resonance amplification after which controllers must be designed if 
conventional methods ar.e followed. -However, resonance amplifications caused 
by small or even zero damping have no effect if a system is reactively stable . 

The proposed method was studied for the Skylab case, where the lag of 
the CMGs was compensated with rate gyros and angular accelerometers such 
that reactive stability was obtainable. However, the case has been simplified; 
e.g., by neglecting sensor lag, CMG cross -coupling, and sampled signal pro- 
cessing. This negligence may be permissible, but generally must be evaluated 
by a more detailed analysis than presented here. 
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Figure 10 . Modular space station configuration. 


The algebraic structure of reactive stabilization leads to a build- 
ing block approach and stability can be predicted by rather simple 
subsystem analyses. The method is ideally suited for syntheses of control 
systems, which is important for multivariable cases with high degrees of 
freedom. The given definitions and theorems are related to the frequency 
response and result in simple criteria similar to Nyquist plots. The method 
promises to simplify design, to untangle different control areas, and to reduce 
overall cost. 

George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Marshall Space Flight Center, Alabama 35812, October 6, 1971 
908-10-08-0000 
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